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Abstract 

We present a numerical and theoretical study of the zero-dispersion limit of the 
focusing Zakharov-Shabat hierarchy, which includes NLS and mKdV flows as its second 
and third members. All the odd flows in the hierachy are shown to preserve real- valued 
data. We establish the zero-dispersion limit of all the nontrivial conserved densities and 
associated fluxes for these odd flows for a large class of real-valued initial data which 
includes all "single hump" initial data. In particular, it is done for the "focusing" 
mKdV flow. The method is based on the Lax-Levermore KdV strategy, but here it is 
carried out in the context of a nonselfadjoint spectral problem. We find that after an 
algebraic transformation the limiting dynamics of the mKdV equation is identical to 
that of the KdV equation. 
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1. Introduction 



1.1. Background. The asymptotic analysis of nonlinear wavepackets, extending the 
Fourier-theoretic treatment of the linear case, has developed into a prominent area of 
activity in both Analysis and Applied Mathematics. The first systematic approach to 
this problem is due to G. Whitham [W], who, in the 1960's, found various mechanisms 
(variational principles, averaged conservation laws) for formally deriving modulation 
equations which describe the evolution of parameters in a family of exact travelling wave 
solutions to a nonlinear equation. This is in analogy to the classical geometric optics 
approach for linear equations. There the asymptotic evolution of a linear wavepacket is 
described in terms of a deformation of parameters in the family of plane wave solutions 
corresponding to the dispersion relation of the equation. 

It was fortuitous that at the same time that Whitham was formulating his theory, 
the subject of completely integrable partial differential equations (PDE's), or soliton 
theory, became widely known. The latter concerns a class of physically significant 
PDE's with large families of exact solutions, including multiphase waves, as opposed to 
just travelling wave subfamilies. This provided a rich and fertile testing ground in which 
to explore the techniques that Whitham had developed. Whitham studied a wide vari- 
ety of these equations himself in the travelling wave setting: Korteweg-deVries (KdV), 
nonlinear Schrodinger (NLS), sine-Gordon (SG), etc. The extension of Whitham anal- 
ysis to the class of multiphase solutions was accomplished by Flaschka, Forest and 
McLaughlin [FFM] initially for KdV. This was still in the setting of wavepacket anal- 
ysis. The first work to place Whitham's theory within the framework of an explicit 
singular limit problem was due to Lax and Levermore [LL2] who investigated the zero- 
dispersion limit of the Korteweg-deVries (KdV) equation. There one studies the limit 
as e — > of the conserved densities for the scaled KdV equation 



The limit is strong and given by the solution of the Hopf equation, which is hyperbolic, 



so long as its solution is classical. After the breaktime the limit is weak due to the 
development of regularizing small wavelength oscillations with an amplitude of order 
unity; thereafter its evolution is no longer governed by the Hopf equation (1.2a). 

In their seminal paper, Gardner, Greene, Kruskal and Miura [GGKM] showed 
that the KdV equation is completely integrable using the inverse scattering transform 



(1.1a) 
(Lib) 



d t u e — 6u e d x u e + e 2 d xxx u e = , 
u e (x,0) = u in (x) . 



(1.2a) 
(1.2b) 



d t u + 6u d x u = , 
u(x,0) = u in (x) , 
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associated with the selfadjoint Schrodinger operator 



(1.3) 



C s = -e 2 d xx + u e . 



Lax and Levermore [LL] realized that this inverse transform would enable one to reduce 
the asymptotic analysis of nonlinear KdV wavepackets to classical short wave asymp- 
totics. They analyzed the limiting behavior of the scattering and inverse scattering 
transform using a WKB analysis of (1.3) and a kind of steepest descent argument to 
obtain a characterization of the (weak) limits in terms of the solution of a variational 
problem. The solution of this variational problem was then constructed through the 
solution of a Riemann-Hilbert problem. Venakides [V] has analyzed the microstructure 
of the limiting solutions, bridging the gap to the local approach of modulation theory 
developed in [FFM]. These results are surveyed in [LLV]. 

Extensions of the Lax-Levermore strategy to other integrable PDE's have been 
pursued, perhaps the most notable being a recent analysis [JLM2] of the defocusing NLS 
hierarchy. In this paper we are going to investigate the asymptotic behavior of a part 
of the hierarchy of commuting flows associated to the focusing nonlinear Schrodinger 
(NLS) equation (1.4). A significant difference between the asymptotic analysis of this 
last equation and that of the KdV or even the defocusing NLS hierarchy is that the 
inverse transform of the former is based on a non-selfadjoint scattering problem (1.5 
+) whereas that of the latter is based on selfadjoint scattering problems. At a formal 
level this has the consequence that the Whitham modulation equations for focusing 
NLS can become elliptic and therefore ill-posed for general initial data. 

This non-selfadjointness also presents an obstacle for knowing how to even begin a 
Lax-Levermore analysis. The work in this paper makes a first step towards confronting 
this obstacle. We do not try to tackle the focusing NLS equation proper; rather, 
we expand our view and look at the entire hierarchy of flows. By abandoning the 
even members of the flow, and the possibilty of complex initial data, we see that the 
remaining members are a hierarchy in their own right, and that their limiting behavior 
is identified with that of the Korteweg-deVries hierarchy in the same limit. Numerical 
simulation suggests that this is a sharp result: after the initial shock, there are 0(1) 
differences in the two flows which go weakly to zero as the dispersion is taken to zero. 

1.2. The NLS and mKdV Equations. The one dimensional cubic nonlinear 
Schrodinger (NLS) focusing (+) and defocusing (— ) equations 




with initial data in the amplitude-phase form 



(1.4b) 
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were first solved by Zakharov and Shabat [ZS] for some natural far field boundary con- 
ditions of the initial data A(x) (positive) and S(x) (real) by using the inverse scattering 
transform associated with a Dirac operator (ZS operator) [ZS1, ZS2]: 



(1.5) C = 



—ie d x ±iu e 
iu e ie drr 



This ZS operator is selfadjoint if in the defocusing case, but nonselfadjoint if in the 
focusing case. The parameter e is introduced into the evolution and scattering problems 
with the zero-dispersion limit scaling; namely, replacing d x and dt by ed x and edf. The 
choice of phase scaling in the initial data is the natual balance for obtaining a non- 
trival momentum density in the limit. This scaling, both in the data and the evolution 
equations, is completely analogous to the linear Schrodinger case of quantum mechanics, 
where the e is taken as a nondimensional version of Planck's constant, h. In the linear 
case, there is no specific scattering operator to contend with. 

For the nonlinear case, the solution of the initial-value problem centers around the 
eigenvalue problem 

7(D 



£f = Cf, where / = 



where now both t and e play roles as parameters, and are not explicitly differentiated. 
For convenience, we may sometimes represent the eigenvalue ( by its real part £ and 
imaginary part rj. 

The NLS equations are the second members, and first interesting members, of a 
sequence of flows associated with £ (1-5) [FNR]. The third members are the modified 
KdV equations (mKdV), 

3 e 2 

(1.6a) d t u e ± -u c2 d x u e + —d xxx u e = , 

or, more generally, their complex forms 

3 e 2 
(1.6b) d t u e ± -\u e \ 2 d x u e + -d xxx u e = 0, 

As members of the focusing (+) and defocusing (— ) NLS hierarchies, they share the 
same respective ZS operator and the same respective conserved quantities. 

There are in fact an infinite number of conserved densities for the NLS hierarchy. 
The first two of them are mass density p e and momentum density [f defined by 

(1.7) p £ = K| 2 , if = -i e -{u"d x u" - u"d x u") . 
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The corresponding conservation laws of the focusing (— ) and defocusing (+) NLS equa- 
tions are 



d t p e + d x p e = 0. 



(1.8) 



d t n e + d x 



A* 1 e2 



-jd x [p e d xx log(p e )] ; 



while those of the focusing (+) and defocusing (— ) mKdV equations are 

d xx p 



(1.9) 



where 



d tl f + d x ( ± ? P v 3 ^ 



\ if 1 I 1 ° r 



4 p e 
3 



„ e d x p £ d xx p e (d x p e ) 2 

TV = -—d x p —p + — — — p . 

p e p e 2p e 



The solutions of either the NLS equations or mKdV equations can then be deter- 
mined by p e and p e up to a constant phase. The other conserved densities can be also 
represented in terms of p e and p e ; for example, the energy density e e for the focusing 
(— ) and defocusing (+) NLS hierarchy, defined by 

(1-10) e e = ^W| 2 t\\u"\\ 

can be determined if p e and p e are known because 

e 1 \p e \ 2 e 2 \d x p e \ 2 1 e2 
2 p £ + 8 p- + 2 P ' 

There are recursion formulas (see equation (7.5)) from which all of the NLS conserved 
densities can be explicitly generated by iteration. 



1.3. The Zero-Dispersion Limit. We want to understand the behavior of solutions 
of the ZS hierarchy when e — > 0. Formally, a natural guess is that the limits of p e and 
p e , denoted by p and p, would satisfy the following equations associated to the focusing 
(-) and defocusing (+) NLS (1.4a) 



(L11) ft/W^T^O; 



d t p + d x p = 0, 
i 2 1 
p 2 
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or the following equations for focusing (+) and defocusing (— ) mKdV, (1.6), 



(1.12) 



d t /i + d. 



dtp + d. 



X 



x 




3^ 

4 p 
4^2 




with initial values which are the limits of the initial p £ and p £ . We call them zero- 
dispersion limits because the formally small terms are dispersive. We call them semi- 
classical limits because these limits in general exist only in a weak sense. 

Consider the equations for p and p given by (1.12) with the focusing sign. They 
are the formal limits of the first two conservation laws specialized to the focusing mKdV 
flow (1.9). In the case that £t = 0, the equation for p is nothing but the Hopf equation 
(1.2a) which governs the zero-dispersion limit of the KdV equation (1.1a) near the 
initial time, except for a trivial scaling: the time variable t in (1.12) is replaced by At 
in (1.2a). 

With a brief calculation, one can find that the systems (1.11) and (1.12) are not 
always well-posed for initial value problems. They are both hyperbolic for the defocus- 
ing case, which is well-posed; but both elliptic for the focusing case, which is ill-posed 
for general initial data. This has been a principal obstacle to the further study of 
the dispersive limit of the focusing hierarchy. However, as we shall describe, there is 
a subset of the flows, namely the odd flows such as (1.12), which will propagate real 
data (S = 0) preserving its reality. (This invariance is decidedly untrue for the general 
even flow; for instance, in the focusing NLS flow itself.) We will demonstrate how the 
analysis of the dispersive limit, on this invariant subspace, can be accomplished. 

Motivated by the results presented here, Kamvissis, McLaughlin, and Miller [KMM] 
very recently undertook a careful analysis of the dispersive limit of focusing NLS for 
a very particular initial condition: A = sech(x) and S = 0. For this data one has ex- 
plicit formulae for the WKB scattering data. Using newly developed Riemann-Hilbert 
problem techniques, they are able to show that the asymptotic behavior of this partic- 
ular evolution problem is correctly described by the Whitham equations, despite their 
ellipticity. Given this development, it now seems probable that this conclusion may 
hold more generally for analytic data. 

1.4. Outline of the Paper. In Section 2, we display several numerical experiments 
to compare the different situations: well-posed and ill-posed system, complex valued 
and real valued system as well as the dynamics among the KdV, focusing mKdV and 
the defocusing mKdV initial value problems with the small parameter e. In Section 
3, we describe the direct scattering data of the focusing NLS equation given by (1.4) 
and its hierarchy. In Section 4, we describe, via WKB asymptotics, the class of "single 
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lobe" real initial data; this asymptotic description can be given in terms of reflectionless 
(multi-soliton) potentials. We analyze the formulae for iV-soliton potentials and deduce 
asymptotic upper and lower bounds on these expressions as well convexity properties 
that enable us to establish a KdV-like limit. In Section 5 we present the rigorous deriva- 
tion of this limit and its relation to a fundamental maximization problem. We then 
show how the densities and fluxes for the odd flows on real initial data can be directly 
expressed in terms of the maximizer of this problem. Section 6 then describes the lim- 
iting dynamics of these flows in terms of a Riemann-Hilbert boundary value problem 
for the maximizer. Finally, in Section 7, we introduce a "loop algebra" description of 
the phase space for the NLS system which emphasizes the Lax pair structure of this 
system. In this setting it becomes easy to understand the behavior of the odd flows 
when one restricts to real initial data. One can identify an invariant subspace for these 
odd flows, even at finite e, which further collapses, formally as e — > onto a space 
isomorphic to the KdV dispersionless densities. This gives a different approach to the 
conclusions found in the previous sections. Finally, we make some concluding remarks 
in Section 8. 



2. Numerical Experiments 

Having introduced the problem, we now present a few numerical calculations that 
will serve to illustrate and motivate the forthcoming analysis. We will be considering 
the real initial value problem for focusing mKdV. 

Figure (2.1) presents a sequence of three snapshots of the evolution of the focusing 
mKdV flow (1.6a-b), for the real, periodic initial data 

(2.1) u(x,0) = cos 2 (x) , 

with e = 0.04. In these figures one can see quite clearly the formation of a region of 
oscillations, of bounded amplitude, within the profile. A fourth panel in this figure 
shows a space-time contour plot of the evolution. For small times one sees that this 
evolution appears to move along straight line characteristics and that the oscillations 
begin to set in where characteristics focus and begin to form an envelope. Indeed this 
evolution very closely matches that of the invsicid Burgers equation (1.6a with e = 0) 
for the same initial data. From the method of characteristics one calculates the break 
time for this dispersionless equation to be ts ~ 0.49, which agrees well with where 
one sees the onset of oscillations. After this time, the solution develops a modulated 
periodic wave train (a so-called "Dispersive Shock") which propagates away from the 
initial kink. Eventually the shock regions collide creating more complex modulational 
structures. 

Figure 2.2 shows the evolution of the same initial data as e is varied between 
e = 0.08, e = 0.04 and e = 0.02. Notice the structure of the shocks as e is varied. 
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It appears that the width, height, and speed of the modulational envelopes remain 
essentially the same as e — > 0, while the oscillations within a given envelope are roughly 
proportional to 1/e. Because of this increasingly oscillatory structure, one does not 
expect to see pointwise convergence, as e — > 0, of the conserved densities (1.7, 1.10) 
mentioned in the previous section. 

However, Figure (2.3) indicates that local averages of these densities might have a 
regular limit: the anti-derivatives of the mass density p e and of the energy density e e 
appear to be converging pointwise to a limit with smaller e producing higher frequency 
but smaller magnitude deviations. Note that fj, e (x,t) = e 1m (u*u x ) = for real initial 
data in the odd flows. 

We conclude this section by mentioning how these simulations were done. Time 
was discretized using implicit midpoint or implicit-Runge-Kutta based on the fourth 
order Gauss points. Space was discretized pseudo-spectrally, with the nonlinear terms 
computed without de-aliasing, which was not necessary because we had adequately 
resolved the spectrum of the nonlinear terms. All these calculations were done in 16- 
digit precision, and the simulations preserved the first three conserved quantities to 
7-digit accuracy. 
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DRAFT 



u(x) for t=0.5 



DRAFT 



u(x) for t=2.0 



DRAFT 



u(x) for t=1.0 



DRAFT 



x-t-u contour plot 



Fig 2.1 Time-slices, u(x, t)\t=o. 5, 1.0,2.0, and contour plot (lower right) u(x, t) = const., 
< x < 27r,0 < t < 4. u(x,0) = cos 2 x. The horizontal axis is x, and the vertical axis 
is either u (slice), or t (contour plot). 
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DRAFT 



DRAFT 



e = 0.08, t = 1.5 



e = 0.08 contour plot 



DRAFT 



DRAFT 



e = 0.04, t = 1.5 



0.04 contour plot 



DRAFT 



DRAFT 



e = 0.02, t = 1.5 



0.02 contour plot 



Fig 2.2. Contour plot, u e (x,t) = const., < t < 1.5 (left). Slice, tt e (x, t) \t=i.5, 
(right). Here, e = 0.08 (top), 0.04 (center), and 0.02 (bottom), and u e {x, 0) = cos 2 a;. 
The horizontal axis is x, and the vertical axis is either t (contour plot) or u (slice). The 
graphs have the same scales for comparison^ 



DRAFT 



d- 1 [M 2 ] 



DRAFT 



\u\ 



Fig 2.3 Primitives of the first (left) and third (right) conserved densities. The cases for 
for e = 0.08 and e = 0.02, are overlaid, with the wider oscillations corresponding to large 
e. Here, the pseudo-inverse d~ x is a periodic normalization of J : each discrete Fourier 
mode of the density is multiplied by jr, except for mode 0, which is left unchanged. 
As e decreases, both the wavelength and the magnitude of the deviating oscillations 
decrease. This suggests that the limit for the primitives will be classical. 



3. The Focusing Zakharov-Shabat Hierarchy 



3.1. The Scattering Transform. Zakharov and Shabat [ZSl] solved the Cauchy 
problem for the focusing NLS equation (1.4) with initial data A(x) and d x S(x) that 
decay sufficiently rapidly as \x\ — > oo. They used the inverse scattering transform 
associated with the nonselfadjoint Dirac operator 

—ied x iu \ 
iu ied x J 

Here we collect the relevant facts regarding their theory. Throughout this section e > 
will be considered to be fixed, and hence, no explicit e dependence will be indicated. 

The solution strategy centers on the spectral problem 
(3.1b) Cf = U where f = ( f L 



(3.1a) 
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Given u(x,0), the asymptotics of the eigenfunctions f((,x, 0) as \x\ — > oo, referred to 
as the scattering data, can be calculated in principle. The evolution of the scattering 
data is governed by simple formulas and u(x, t) is then obtained from the knowledge 
of the large |x| asymptotics of /(£, x, t) using inverse scattering theory. 

More specifically, if u(x,0) has the form (1.4b) where A(x) and d x S(x) decay 
sufficiently rapidly as x — > ±00 (say, faster than any power of x), then the L 2 spectrum 
of £ consists of the whole real axis, comprising the continuous spectrum, along with a 
finite set (possibly empty) of simple eigenvalues in the complex plane. It is easily seen 
that if / and ( satisfy (3.1) then so do / and ( where 



(3.2) 



/ 



/(2) 
_/(!) 



We will therefore only consider (1, • • • , ( N: the N eigenvalues of C located in the upper 
half of the complex plane. 

To describe the asymptotic behavior of the L 2 scattering data we first introduce 
the Jost fundamental solution F(( 7 x). For each ( G C, this matrix solution is uniquely 
defined by the following two conditions: 



i) F(C,x)exp(-zo- 3 Cx/e) 



I 



x 



+00 



where 03 



ii) 



-1 



sup \\F((, x)\\ < 00 . 
xeR 



. This Jost solution is meromorphic in the open upper half-plane, 



as well as in the open lower half-plane, and has limiting values, 

F±(X, x) = lim F(X ± ze, x) , 

from above and below respectively, on the real axis A G R, which are continuous in A. 
These two limiting values are related by a "jump matrix" as follows: 



F + (A) = F_(A) 



l + I^A)! 2 
R{\) 



R{\) ^ 



where the function -R(A) appearing in the jump matrix is called the reflection coefficient. 
It is complex- valued and depends only on A. 

The first and second columns of F((, x), respectively denoted by fi and f^, have 
the following asymptotics in (: 



fi((,x) = exp( 
/ 2 (C, x) = exp 



\ e 
-i(x 



+c(r i ) 

+ o(c 1 ) 



for C - 
forC 



00 , 



00 , 
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The eigenvalues Q in the upper half-plane are determined by the condition that these 
two solutions should be linearly dependent. (These are also precisely the locations 
where the Jost function in the upper half plane has a pole. Away from these points 
the Jost function is analytic.) 

At an eigenvalue, the proportionality between these two solutions 

fi(Q,x) =exp^- ^Jf 2 ((j,x) 

is specified in terms of Xji which is referred to as the norming exponent. 

The inverse theory prescribes that the fundamental scattering data consists of the 
reflection coefficient R(X), the N eigenvalues Q and the N norming exponents Xj- 

As u(x,t) evolves according to the focusing NLS equation (1.4), the eigenvalues 
{Q} remain unchanged while the time dependence of the other scattering data is 

(3.3) Xj (t) = Xj + K 2 jt , R{\ t) = R{\) exp 

Given R(X), Q, and Xj computed from the initial data u(x, 0), the solution u(x,t), as 
well as all of the conserved densities defined in the next subsection, of the NLS equation 

(1.4) , are then determined by inverse scattering from R(X, t), Q, and Xjif) given by 
(3.3) [ZS1, KMM]. 



2i\H 



3.2. The Hierarchy. The complete integrability of this cubic Schrodinger equation 
implies the existence of an infinite family of independent, conserved quantities [ZS1], 

/oo 
p m dx for m = 0, 1, • • • , 
-oo 

which are in involution with respect to the Poisson bracket, 

n rrr rr X 1 f°° fSH m SH n 8H m 5H n \ 

= {H rn ,H n } = — — — — —— \dx. 

ie j-oo \ ou ou ou ou J 

The first three densities p m are essentially those given by (1.7) and (1.10): 

|2 



(3.4) 



Po = ~P = -\u\ 
pi = p = — i- (ud x u — ud x u) , 
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A recursion formula that generates all the conserved densities is given by (7.5). 

All of the H m are Hamiltonians that generate flows which commute with the 
focusing cubic NLS flow (1.4a), the so-called focusing Zakharov-Shabat (ZS) hierarchy. 
Letting t m denote the time variable associated with the m th flow, its evolution is then 
given by 

ied t u = - m for m = 0, 1, • • • . 
8u 

The to how is just a phase rotation, the t\ flow is just a spatial translation, the flow 
is given by the focusing NLS equation (1.4a) with t = t 2 , while the £3 flow is that of 
the complex mKdV equation 



d t3 u + -\u\ 2 d x u + e 2 ^d xxx u = d 



Every H n is conserved by each flow; their densities satisfy the local conservation laws 

to r\ dtoPn-1 =0, 

(3.5) for m,n = 1,2, ••• . 

Ot m Pn-l + Ox^m,n = U , 

Here n m ,n is the flux for the (n — l) st conserved quantity under the m th flow. 

Because all these flows commute, they may be solved simultaneously for u(x,t), 
having the form of the initial condition (1.4b), where t = (to, t±, • • • ) such that all but 
finitely many t m are zero. Associated with each t is a polynomial p( ■ , t) defined by 

00 

(3.6) p((,t)=J2tm( m . 

m=0 

The simultaneous evolution of the scattering data is then given by 

(3.7) Xj (t) = Xj+ip{M , A(A,t) = tf(A)exp(^M) , 

and u(x,t) is determined by inverse scattering. 

In Section 7 we show that the odd flows in this hierarchy preserve real symmetry. 
More specifically, we show that u(x : t) is real whenever w(x,0) is real and t m = for 
every even m. We show moreover that when u is real that p n = for every odd n. 

The scope of the zero-dispersion limit for the NLS equation can then be enlarged to 
consider the solution u e (x,i) of the whole hierarchy that satisfies the initial condition 

u e (x, 0) = A(x) exp^^- 1 
14 



where A(x) is positive, S(x) is real, both are independent of e, and A(x) and d x S(x) are 
Schwartz class. The goal is then to determine the limiting behavior of all the conserved 
densities p e n and fluxes \f m n associated with the entire focusing ZS hierarchy of flows 
as e tends to zero. In this paper we achieve this goal only for the case when u is real 
and we restrict to the odd flows of the hierarchy. 



4. Analysis of the Scattering Transform 



4.1. Asymptotic Analysis of the Initial Scattering Data. For initial data 



u e (x,0) = A(»exp -S(x) , 



% 



e 



of the class described at the end of the previous section, one can asymptotically deter- 
mine the initial scattering data by a WKB analysis. 

The L 2 eigenvectors as well as the reflection coefficient, -R(A), for the Zakharov- 
Shabat eigenvalue problem (3.1) were defined in Section 3. The WKB ansatz is that 
the eigenvectors have a leading order asymptotic (in e) behavior of the form 

r a> = r c» e J M + . . . , r m = r m ew (>^±M) + . . . , 



where A 2 \ and 9 do not depend on e. (There is a similar ansatz for the Jost 
functions.) When this ansatz is substituted into Cf e = (f e , the leading order gives 

( C + {d x 6 - \d x S) iA WrW\ 

V • A t - {d x 6 + \d x S) ) \rV> ) U - 

This implies that 

V iA C - {d x e + \d x s) ) 

= (C - \3 X S) 2 - (d x 9) 2 + A 2 . 

Suppose now that the initial value of u e is real. That means u e (x, 0) = A(x) and 

(dj) 2 = (A 2 + C 2 ) . 

If the WKB asymptotics are valid then, in order for f e to be an eigenvector, d x 9(x) 
should be pure imaginary for sufficiently large values of \x\. It follows from the preceding 
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equation that ( should then be pure imaginary. Thus, we would expect the point 
spectrum to converge, as e — > 0, onto the imaginary axis. 



In fact there is a natural class of initial data, the so-called "single lobe" potentials 
for which the eigenvalues are always pure imaginary independent of the value of e. By a 
single lobe potential we shall mean a real- valued potential for which A(x) is nonegative, 
smooth and integrable on R and moreover A(x) should be non-decreasing on the left of 
x = and non-increasing on the right. It has recently been demonstrated [KS] that for 
such data the eigenvalues are always pure imaginary. In fact these authors show that 
this result holds under a somewhat weakened definition of single lobe, but the above 
will suffice for our purposes. For such data, A(x) will have two turning points, x^(r)) 
(Fig. 4.1). We will assume that A(x) has maximum height 1 . In this setting the WKB 
ansatz implies that the asymptotic scattering data should satisfy: 

(1) The reflection coefficient -R(A) = 0. 

(2) The eigenvalues are independent of time. They are purely imaginary with Q = irjj. 
The density of states, for the positive imaginary eigenvalues, is 



1 d r i d r x +(v) 

^fa) = - :T Re / V A2 (y)-V 2 dy = --r- \ y / A 2 (y)-r 1 2 dy; 

the eigenvalues are then defined by 



flmax 

Jrij 



for j = 1, 2, • • • , iV e and t]i < r)i < ■ ■ ■ < t]n^ where 

(3) iV e is the total number of eigenvalues in the upper half plane and is determined 
by 



iV e = 



(4) The initial norming exponents can be approximated by 

Xj(0) = Xj = Vj x + (Vj) + [ (vj ~ Jv]-A 2 (x)) dx 
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DRAFT 



Fig 4.1 The initial data A(x). Note its critical value r\ max and the indicated defining 
relations for x^ 1 = x ± (rj). 

Thus, if one just considers the odd flows, i.e. t = (— 1\, £3, . . . (— l) fc t2fc-i, • • • ) then 
all the scattering data are real. The formula for Xj given by (3.7) can be rewritten as 



We observe that this asymptotic data is the same as that for the Schrodinger operator 
(1.3) with A 2 (x) playing the role of Ui n (x) [LLV]. Henceforth, we consider the Cauchy 
problem for focusing NLS (1.4a) with initial data corresponding to the above asymptotic 
scattering data. 

4.2. Reflectionless Approximation. Motivated by the previous section, we choose 
to replace the exact real- valued initial data A by the real- valued reflectionless potential 
A e corresponding to the WKB scattering data given above. This allows us to use 
the reconstruction formulas for reflectionless potentials [FNR] that give the associated 
conserved densities and fluxes by 



where the so-called tau-function r e (x,t) is given by the A^ e x A^ e determinant 



oc 



(4.2) 




(4.3) 



Pn-i(x, t) = e 2 d xtn log(r £ (x, t)) , 
/C,n(M) = -e 2 d tm t n log(r £ (x,t)) , 



(4.4) 



r £ (x,t) = det (I + e 2 G e (x,t)G e (x,t 



))• 
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Here I is the identity matrix while the N £ x N £ matrix G e has the form 



j,k=l 



(4.5) G*(x,t) = ( — U= e ^(«i,^) + <^*X 

. - Cfc V 

where 

(4.6) a(r],x,t) = -rjx + Xj(t) ■ 

All the eigenvalues {Cj}i he on the positive imaginary axis and are bounded as 
follows: 

(j = iVj, < r){ < rj e 2 < ■ ■ ■ < ?]% < 1 , 
then E is real for all the odd flows (t = ti, £3, • • • ) and 

% 1 



- Cfc + Vl 
From this, we see that G e is real, positive definite and symmetric. 
In this setting we define 



B(x,t) = G e (x,t) = G e (x,t). 



4.3. Convexity of the Tau-Function. We will make use of the following property 
of the tau-function. 

Lemma 4.1. For each e > the function 

(x, t) 1— > t(x, t) is convex . 

Proof. Set W = e 2 log det(I + e 2 B 2 ). It is sufficient to show that W is convex in x for 
general real parameters rjj and Xj- This is true because convexity is a property along 
line segments, and examining W along a general line, parameterized by x, of the form 
P(x) = (x(x), t(x)) = Po(l — x) + P\x results in the function 

W(x) = e 2 log det(7 + e 2 B(x) 2 ) . 

The convexity of W(x, t) follows from the convexity of W(x). For the remainder of the 
argument, we drop the tildes and denote differentiation in x with a prime ('). 
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Using the definition of 5, we find B' = — j(AB + -BA), where A is the diagonal 
matrix of the rjj. From this and the general facts that det(a)' = det(a)( tr(a _1 a') and 
tr(a&) = tr(6a), we find by direct calculation that 

(4.7) W = -4e tr((J + e 2 B 2 )~ 1 e 2 B 2 A S j . 

It is more convenient to write W as 

W' = -4e tv^I+(e 2 B 2 y i y 1 A ) 
before differentiating to obtain 

(4.8) W" = 4 tr [(A)(AB 2 AA)} + 8 tr [(AB)(ABAA)} + 4 tr [(AB 2 )(AAA)} , 
where 

A=(l+(eB)- 2 )~ 1 (eB)- 2 . 

Both A and B are symmetric positive definite (SPD), and A commutes with B, so 
products of A with powers of B are also SPD. Thus the terms of W" in (4.8) are traces 
of the products of two SPD matrices, each grouped in parentheses in the equation. 

Traces of products of SPD matrices are positive, and so the sum is positive. Thus 
W" > 0, and so W(x) is convex, and, so, W(x,i) is convex, which proves Lemma 4.1. 

We remark that (4.7) for the full system yields 

e 2 a x logr £ (+oo,t) = 0, e 2 d x log r e (-oo, t) = -4e tr(A) . 
One can similarly find e 2 d tm logr e (±oo,t) for any positive, odd number m. 



5. Establishing the Limit 



5.1. The Limit of the Tau-Functions. In this subsection we establish the 
limit, as e tends to zero, of e 2 log(r e (a;, t)) in the topology of uniform convergence 
over compact subsets of R x T Ql characterizing the limit in terms the solution of a 
maximization problem. Here, T Q denotes the subspace of vectors of the form t = 
(— 1±, ^3, . . . ( — l) k t 2 k-i, • • • ) which are real and have finite support. 



19 



We adopt the following notation. Given any N e x N e matrix A indexed by the set 
J e = {1, . . . , iV e }, and given any two nonempty subsets S C J e and T C J e , let A$t 
denote the \S\ x \T\ submatrix of A obtained by retaining only those entries whose row 
index belongs to S and column index belongs to T. We abbreviate Ass as As- Finally, 
for every £ C J e , we let 5" denote the compliment of S in J e . 

The key new result is the following. 
Theorem 5.1. Let 

(5.1) q £ (x,t) = max{e 2 log(det (G e s (x,t))) :ScJ £ }, 

where the maximum ranges over all index sets, with the understanding that 
det(Gg) = 1 when S is the empty set. Then 

(5.2) lim e 2 log(r e (x,t)) -q £ (x,t) =0 
uniformly over R x T Q . 

Proof. Below we will derive the crude bounds 

(5.3) e 2 ^exp(^^) < r ^ t ) < 2 2N * exp ^^i^l^ . 

Given these, by taking their logarithm, multiplying the result by e 2 , and subtracting 
q e (x,t), one arrives at 

(5.4) e 2 2N e log(e) < e 2 log(r e (x, t)) - q £ (x, t) < e 2 2A^ e log(2) . 

Because N e = 0(1/ e) by the WKB ansatz, the limit (5.2) follows. Hence, the theorem 
follows once the bounds (5.3) are established. 

To derive the upper bound on r e in (5.3), first use the fact that G e is Hermitian 
positive to obtain 

r e = det(l + e 2 G e2 ) 

(5.5) V ! 2 

< det(/ + 2eG e + e 2 G e2 ) = [det(l + eG e )] . 

The last determinant in (5.5) is just the characteristic polynomial of eG e evaluated at 
—1 and can be expanded in terms of the determinants of its minor matrices eG s as 

(5.6) det(/ + eG e ) = ^ e' 5 ' det(G £ 5 ) . 

ScJ e 
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Because every principal minor of G e has the same form as G e , each G e s is Hermitian 
positive, and hence, has a positive determinant. The sum in (5.6) therefore contains 
2 N positive terms, each of which is bounded above (for e < 1) by 

(5-7) c ISI de t(G|)<exp(*^), 

The upper bound on r e in (5.3) therefore follows by combining (5.5), (5.6), and (5.7). 
To derive the lower bound on r e in (5.3), first use the fact that 

(5.8) r e = det(J + e 2 B) = det(e 2 B s ) . 

ScJ* 

Because G e is Hermitian positive, so is B = G e2 . Hence, each term of the sum in (5.8) 
is positive. We then use the fact that 

and the fact that G e ss G e ss and G e ss ,G e s , s are Hermitian nonnegative matrices to bound 
each term below by 

det(e 2 S s ) >det(e 2 G^ 5 G £ ss ) 

(5.9) =e 2 l s l[det(G^)] 2 

>e 2 Hdet(^)] 2 . 

Because each term of the sum in (5.8) is positive, the sum may be bounded below by 
the largest lower bound on any of its terms. By (5.1) that is achieved for the S where 

[det(^)] 2 =exp(^)), 

whereby the desired lower bound on r e in (5.3) is obtained, and the proof is complete. 
5.3. Properties of the Maximization Problem. 

Before we can characterize the semiclassical limit of the conserved densities and 
fluxes, more basic facts concerning the maximization problem in (5.1) must be estab- 
lished, including the uniqueness of the density at which the maximum is attained. By 
introducing the scalar product 

1 f 1 

(5.10) {a\(3) = - a{r,)P{ri)dr,, 

n Jo 
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and the integral operator 



(5.11) 



dv . 



The underlying functional being maximized in (5.1) can be recast more abstractly as 
the quadratic form 



where a is again given by (4.6) but here we suppress the rj dependence. The properties 
of these objects that we will need are presented in the following three lemmas. 

The maximization problem (5.1) in this form is then posed over the set A of 
admissible L 1 densities, 



which naturally inherits a topology from the weak-* topology of measures. However, 
the densities in A, being all bounded above by p(ry), are equi-integrable and, hence, 
comprise a relatively compact set in the weak topology of L 1 ([0, 1]). This means that 
the weak-* topology of measures and the weak topology of L 1 coincide on A. We will 
always consider A equipped with this topology. The key fact we will need is stated in 
the first lemma, the proof of which is omitted. 

Lemma 5.2 (^-compactness). The set A is compact. 

The operator L defined by (5.11) and the functional Q defined by (5.12) are well 
behaved over the set A. Indeed, given that ifj is in A, it can be shown [LL2] that Lif; is 
in the class of continuous, odd, functions of ?y, which we denote as C oc id(R)- Moreover, 
it is clear from (4.6) that for every (x,t) 6RxT the function rj i— > a(?y,a;,t) is also 
in C dd(R). Consequently, Q defined by (5.12) takes values in R. Moreover, L and Q 
possess continuity properties that are stated in the next lemma. 

Lemma 5.3 (^-continuity). 

(a) The operator L : A — > C G dd(R) is continuous. 

(b) The functional Q:ylxRxT ^Ris continuous. 

Proof. The hardest part is the proof of (a) which resembles that of Theorem 3.4 in 
[LL2]. The proof of (b) follows from (a) and the continuity of a. // 



(5.12) 



Q(ip; x,t) = §(V | LVO + (a(x, t)\if>). 



^ = {^GL 1 [(0,l)]:0<^(r7))<^)}, 
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Remark. By combining Lemma 5.2 with Lemma 5.3, it is seen that the image of A 
under L is a compact subset of C oc Ld(R.) and that for every (x,t) G R x T Q the map 
ip i — y Q(i/j;x,t) has a compact range and, hence, attains a maximum over A. The 
existence of the maximum for (5.1), which has already been established through the 
limiting procedure of Theorem 5.1, is thereby re-established intrinsically. 

The uniqueness of the maximum in (5.1) follows from the third lemma. 

Lemma 5.4 (strict concavity). 

(a) The quadratic form (ip | Lip) is negative definite over A A. 

(b) For every (x, t) G RxT„, the map ip i— > Q(ip; x, t) is strictly concave over 
A. 

Proof. The hardest part is again the proof of (a) which here resembles that of Theorem 
3.7 in [LL2]. The proof of (b) follows directly from (a). // 

With the above pieces in place, the main result of this subsection can now be 
persented. 

Theorem 5.5 (uniqueness and regularity). For each (x,t) G R x R°° 
there exists a unique ip*(x,t) in the admissible class A such that 

(5.13) q(x,t) = Q(iP*(x,t);x,t). 

Moreover, this maximizing density satisfies the following. 

(a) The map (rr,t) h- > ip*(x,t) is continuous from RxT„ into A equipped 
with the weak L 1 topology. 

(b) The map (x,t) \— > Lip*(x,t) is continuous from R x T Q into C G dd(§) 
equipped with the uniform topology. 

(c) The map (x,t) i— > q(x,t) is differentiable from RxT„ into R with its 
continuous partial derivatives given by 



(5.14a) d x q(x,t) = -(ri\il>*(x,t)), 

(5.14b) d t2n _ 1 q(x,t) = (v 2n - 1 \r(x,t)). 

Proof. Lemma 5.4 states that for every (x, t) the functional Q(-; x, t) is strictly concave 
over A, thereby insuring the uniqueness of the point in A at which the maximum in 
(5.1) is attained. Denote this point as ip*(x,t). 
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To prove the continuity asserted in (a), let (x,t) be any point and {(i^tj} any 
sequence converging to that point in R x T Q . Because, by Lemma 5.2, A is compact, 
the sequence {tfj*(x k , t^)} has cluster points, all of which lie in A. Let -0* denote one 
such cluster piont. Upon passing to a subsequence if necessay, part (b) of Lemma 5.3 
implies 

(5.15) lim q(x k ,tk) = lim Q(ip*(x k , t k ); x k , t k ) = Q(ip*; x, t) . 

On the other hand, the continuity of q implies 

(5.16) lim q(x k ,t k ) = q(x,t) . 

k^oo 

Comparing (5.15) and (5.16) shows that Q(ip*;x,t) = q(x,t), whereby we conclude 
that ip* = ip*(x,t). Hence, the original sequence {i()*(x k , t k )}, having ip*(x,t) as the 
only cluster point, must converge to tp*(x,t). The continuity asserted in (a) follows 
immediately, while that of (b) does so after invoking the part (a) of Lemma 5.3. 

Now turn to the differentiability (c). It suffices to establish (5.14b). Let t and t' 
differ only in the coordinate tm-\- For every rj in A, a direct calculation then yields 

(5.17) W; x, t') - W; x, t) = {r, 2 ^ 1 \ V ) (^ n _ x - t 2n _ x ) . 

However, because tjj* maximizes Q, one has the general two-sided inequality 

(5 18) Q(^( x ^^ x ^') - W*(a:,t);x,t) < q(x,t') - q(x,t) 

< Q(ip*(x,t');x,t') - Q(tp*(x,t');x,t) , 

which, when combined with (5.17), gives 

1Q1 {V 271 - 1 | P{x, t)) (4 n -i " t 2 n-i) < q(x, t') - q(x, t) 

^(t/ 2 "" 1 1^(^,0) (4.-1 -*2„-l). 

The result now follow by the continuity of rj*. / / 



It will prove useful that the solution tfj*(x,t) of the maximization problem can be 
characterized in terms of variational conditions. 

Theorem 5.6 (variational conditions). If rj E A then rj = rf(x,t) if and 
only if rj satisfies the variational conditions 



(5.20) V = 



where a(rj, x, t) + Lifj < , 
4>{rj) where a(rj, x, t) + Lifj > . 
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Proof. The proof follows that of Theorem 3.12 in [LL2]. It uses the continuity of 
a + Lip in 77 to assert that the conditionals in (5.20) define open sets. / / 

5.4. The Limit of the Densities and Fluxes. By combining the differentiability 
result of Theorem 5.5 with the following elementary, but nontrivial, lemma will yield a 
strengthening in the sense of the convergence for the tau- functions in (5.2). 

Lemma 5.7 (converging derivatives). Let {h n } be a sequence of differ- 
entiable convex functions over R°° such that h n (x) — > h(x) uniformly over 
compact subsets of x G R°° and h is differentiable. Let d s = x • <9 X denote the 
directional derivative of x in a given direction x G R°°. Then 

(5.21) cUn(x) -+d a h(x), 
uniformly over compact subsets of x G R°°. 

In particular, we can now give the main result of this section. 

Theorem 5.8 (limit of densities and fluxes). The limit of the tau- 
function r e is given by 

(5.22) hme 2 ^ logr £ (x,t) = d t2n _ iq (x,t) = ( V 2n ~ l | t)) , 

uniformly over compact subsets of (x, t) G RxT„. The densities p e n and fluxes 
H e myn have the distributional limits 

(5.23a) V'{dx)-\impl n _ 2 = ^(r? 2 ^ 1 | V*) , 

(5.23b) V(dt2m-i)- lim /4m-i,2»-i = -ft am -i (^ 2n_1 I P) ■ 

Proof. Assertion (5.22) follows directly from (5.2) of Theorem 5.1 (tau-function limit) 
and Theorem 5.14 (uniqueness and regularity) upon applying Lemma 5.7 (converging 
derivatives). Assertion (5.23) then follows by the usual integration by parts argument. 

// 

Remark. The significance of (5.23) is that one can now pass to the limit in the local 
conservation laws 

(5-24) d tam _ 1 P2 n -2 + dxf4m-l,2n-l = Q > for TO, 71 = 1, 2, • • • . 



6. The Limiting Dynamics 
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By the results of the previous section we have seen that the limit related to r can be 
converted to a variational problem. This can in turn be converted to a Riemann-Hilbert 
boundary value problem which can be solved with the use of the Hilbert transform. 
We will list some related results but omit most proofs. The interested reader can find 
more details in [LL] and [LLV]. 

The variational conditions (5.20) for the maximization problem can be expressed 

Lip (if) + a(rj, x, t) < , when ip(rf, x, t) = 0, 
(6.1) Lip(r]) + a(?7, x, t) = , when < ip(r], x, t) < 

Lip(rj) + a(rj, x, t) > , when ip(rj, x, t) = 4>(r)). 

Let / be the interior of the set of (77, x, t) in which equality holds and I be its closure. 
Directly differentiating (6.1) while using formula (4.6) for the function a leads to 

^ 62 ^ Lip x (r]) = 77 , Lip t (r]) = -dtxim*) » when (77, x, t) e I, 

^(r/O = 0, Vt(??)=0, when (77, x, t) £ /, 

where x(r?; t) is defined as Xj(t) if 77 = 77,-. 

These differentiated variational conditions do not have any explicit dependence 
on (x,t), but rather their dynamics, as well as all their memory of the initial data, 
is contained in the set /. Consequently, these conditions have more general validity 
than the variational conditions (6.1). Indeed, the latter vary (although in a sense 
insignificantly) when the initial data are considered in different classes (such as periodic 
or tending to different limits as x — > ±00), while the differentiated conditions remain 
unchanged. 

Making the ansatz that at fixed (x,t) the set J(x,t), defined as 

l(x,t) = Le[0,l) : (77,x,t)e/l, 



consists of a finite union of disjoint open intervals, one can uniquely determine the 
functions ip x (v) an d iptiv) m terms of the endpoints of these intervals. More precisely, 
I(x, t) is assumed to take the form 

(6.3) /(x, t) = [0, /3i) U /3 3 ) U • • • U {fa g , /3 29 +i) 

for some nonnegative integer (7, where the (3i and g depend on (x, t) and 

< Pi < P2 < ■ ■ • < fog+l < 1 • 

The operator L given by (5.11) plays an important role for seeking the unique 
solution of the maximum problem. 
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Extend ip(r)) to be an odd function over the real axis R that vanishes when \rj\ > 1, 
then the ry-derivative of L is the Hilbert transform: 



This observation allows the differentiated variational conditions (6.2) to be transformed 
into Riemann-Hilbert problems through which they can be solved explicitly for ift x and 
ifj t in terms of hyperelliptic functions involving a radical function R. 

This R has the following form 



where the {/^} are the boundaries of the set I(x 1 1) as given by (6.3). 

Note that I{x,i) is the set over which R(rj) is real valued and g is just the genus 
of the associated Riemann surface. 

The variational problem (5.12) is uniquely solvable, say, by ip*. It follows from 
(5.14) and (5.23), for n=l, that all the limits of nontrivial conserved densities along 
odd flows can be represented by the derivatives of ip* in the sense of distributional 
convergence: 



for m = 1, 2, • • • . 

The integrals above are calculated on /(x,t), the particular domain of r] as given 
by (6.3). These {(3j}, as the boundary points of I(x, t), are then governed by a genuine 
nonlinear hyperbolic system under any odd flow. Different flow makes different time 
evolution of the domain I(x, t) through its boundary {A}, and enforces different time 
evolution of the semiclassical limit. 

The dependence of the Pi on (x, t) is then derived from the compatibility constraint 
dti^x = d x tfjt- This reduces to the first order system of hyperbolic equations in the 
Riemann invariant (diagonal) form 





i=l 



p(x,t) = -(r/lVC); 
2 

= g faints). 
p 2m _ 2 (x,t) = (t/I^.J, 



(6.4) 



dt m Pi + Smiiftl, • • • , /?2g+l) 9 x (3i — , 



for i = 1, • • • , 2g + 1 



for any positive odd number m. 
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In this way, the initial value of ip* can be represented in the form 

{0 for x > x + {rf) , 

i: +iv) (A2(y) % 2) i/ 2 d y fOT > x > x-( v ) 
4>{rf) for x < x~ (77) , 

where x~{rj) < x + (rj) are implicitly defined by 

A(x + ) = A(x~) = rj. 



Consider the semiclassical limit along the ts flow, i.e. the mKdV flow (1.9). Note 
that its formal limit equation (1.12) has no classical solution beyond a finite time, say 
t = tb, because of its hyperbolic nature. If p(x,t) is the solution of equation (1.12) 
(focusing) with the initial value p(x, 0) = A 2 (x), then tb is given by 

t b = [max(3AA a; )]~ 1 . 

X 

For t < tb, i.e. g = 0, we have that f3j(x,t) = p(x,t), an exact solution of (1.2), 

and 

L 2 - lim/9 e (x,t) = 0i(x,t) 

uniformly in t. 

For t > tb, the limits of all conserved densities and their fluxes exist in the sense 
of L 2 weak convergence. They are determined by the 2g + 1 functions {fa} which are 
governed by a strict genuinely nonlinear hyperbolic system (6.4) and will break down 
again after a finite time. 



7. An Invariant Subspace for the Hierarchy 



Loop algebras provide a natural "phase space" in which to study the focusing 
Zakharov-Shabat (ZS) hierarchy. We will see that, when Ui n is restricted to be real, 
the correspondingly restricted phase space is invariant under the action of the odd 
flows of the hierarchy. Reduction to this invariant subspace provides a framework in 
which to geometrically understand the dispersionless limit of these odd flows onto the 
corresponding limit for KdV. 



7.1. Loop Algebra Structure. We begin by summarizing the loop algebra structure 
for these soliton equations. Further details may be found in [FNR]. To facilitate our 
presentation, we introduce the following notation for the Pauli basis of sit (2): 
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With respect to this basis the ZS spectral problem can be expressed as 

(7.1) £-( = -eHd x +a£ + [3F-(, 

where u e (x) = a(x) + i(3(x). Upon multiplying this expression by 7i, we get an equiv- 
alent operator, 

(7.2) ed x -aF + (3£ - Qi. 

Observe that the nondifferential part of this operator is a real linear combination of 
the Pauli matrices, and so is a family of elements of su(2). This observation motivates 
the definition of the loop algebra for the hierarchy, which is defined as the vector space 
of all formal Laurent series, Q, with coefficients in su(2): 

oo 

(7.3) Q= C (h n H + a n f + b n £) . 

n= — oo 

The "algebra" structure is a Lie algebra structure induced by the usual Lie bracket, 
[ • , • ] for su(2): if Q = E^-oc C j Q-v and R = T,T=-oc ( k R-k, then 

oo 

^ = E E C[Q-j,R- k ]. 

m= — oo j-\-k=m 

This loop algebra is the natural setting in which to develop the Lax pair formulation 
of the NLS flows. For example, in the case of the original focusing NLS flow if we let 

denote the operator (7.1), then is an element of the loop algebra (7.3) which is 
linear in (. There is another operator, 



ed t + Q 



(2) 



where 



Q(2) =( 2 Q + (Q 1 + Q 2 

C 2 n + c 



( 7 - 4 ) >2^f0 -u\_( ^l 2 ^ 



i-o x u — \u\ 
2 2 1 1 ■ 

is a quadratic polynomial in the loop algebra (7.3) such that the compatibility condition, 

[d x - Q (1 \ d t - Q^} = Of? - Qi 2 ) + [Q (1 \Q {2) ] = , 
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known as a Lax equation, is equivalent to the requirement that u e should solve (1.4a). 



7.2. The Hierarchy of Densities and Flows. In Section 3 we took an approach 
for constructing the hierarchy of flows which is based on a Hamiltonian formalism 
applied to a generating function for the conserved densities. However, this recursive 
construction can be encoded [FNR] in the loop algebra representation as follows. The 
hierarchy of conserved densities, p m = h m+ 2 is recursively defined by 



(7.5a) h n = i J (ue n — ue n ) dx , 

where 

(7.5b) e n = i/2(ed x e n -i - 2iuh n -\) , 

with ho = — 1 and eo = 0. Setting e n = a n + ib n then determines the coefficients, 
a n , 6 n , h n of a loop element (7.3). 



We now introduce the matrix 
(7.6) Q n = 



6 n ih n 



where e n and h n are the quantities we have just defined. One sees immediately from the 
definition of these quantities that Q n is an element of su(2); i.e., it is skew-Hermitian. 
Using these recursive definitions it is also a straightforward exercise, which we leave to 
the reader, to show that the Q n satisfy a hierarchy of coupled ODE's: 

(7.7) d x Q n -[Q 1 ,Q n ] = [Q ,Q n+1 ]. 
This hierarchy can be encoded in a single "generating" ODE 

(7.8) d x Q-[(Q o + Qi,Q] = 0, 

where Q = Qo + C _1( 3i + C _2( ?2 + • • • is a Laurent element of our loop algebra, as 
is again straightforward to check. The coefficient of C _n on the left-hand side of (7.8) 
is equivalent to the nth ODE of (7.7). Note that the linear polynomial in the bracket 
of (7.8), (Qo + Qi = is the nondifferential part of H(C — Q which is equivalent 
to the original Zakharov-Shabat eigenvalue problem. Moreover, if \1> be a fundamental 
solution of the ZS problem, then [FNR] Q = ^H^- 1 solves (7.8). Thus one sees that 
the loop algebra elements have a realization in terms of "squared eigenfunctions" of 
the ZS eigenvalue problem and, therefore, the densities h n and e n can, in principle, be 
read off from the asymptotic expansion in ( of the Jost function for the ZS problem. 
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7.3. Invariance of Real Symmetry for the Odd Flows. Before investigating 
how the loop algebra structure reflects a restriction to real initial data, we need to 
introduce one other ingredient: a natural inner product on the underlying vector space 
of loop elements. For any two loop algebra elements, Q = Y^JL-oc C J Q-j an< ^ ^ = 
Xlfcl-oo ( k R-k, this symmetric form is defined by 

(7.9) (Q )jR ) = _I t r( Q-jR-k 

\j+k=o 

Henceforth it will be assumed that u(x) is real-valued; i.e. that u(x) = A(x) and 
S(x) = 0. Let's use the recursion formulas (7.5) to calculate the first few densities 
under the reality constraint: 

e = 

(7.10) ei = uh = -u 
e2 = l -d x (-u) 

By using this as a base step, we can prove 

Proposition 7.1. If u e R, then e ev E zR, e a £ R and h d = 0. Here the 
subscripts ev or od respectively stand for any subscript that is even or odd. 

Proof. The proposition is proved by induction with the base step given by (7.10). Now 
assume that the proposition has been established for e2 n -2 and hm-2 which are thus, 
respectively, pure imaginary and real. Because, by (7.5b) e2n-\ = § ^d x em-2 + uh-2n-2 
it follows that e2 n -i is real. From this it follows, using (7.5a), that /i2n-i vanishes 
identically. Using (7.5b) again, we have e2 n = \^Qx^2n-\ £ ^R which completes the 
induction. 

Now let Q( ev ^ and Q^ od \ denote any linear combination of loop coefficients of the 
form Q2j and Q2j-\ respectively. (This includes, as a particular instance, loop algebra 
elements which are even in ( and odd in ( respectively.) Then based on the preceding 
proposition one has the following corollary. 

Corollary 7.2. If u e R, then e R-H R-£ and thus in particular has 
pure imaginary components. Q^ od ^ £ R-JF and therefore has real components. 
Moreover, because 7i,S,J-' are an orthonormal basis of su(2) with respect 
to the inner product introduced earlier, one always has JL with 

respect to this inner product. 




ho = -l 
hx = 

h,2 = i j u(—i/2d x u) — u(i/2d x u) dx = -u 2 . 
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Proof. By the proposition, /12.7-1 is identically zero while e2j-i is real. Using this 
in (7.6) one immediately observes that Qij-\ is in the span of T . Because the h- 
coefficients are always real, hi{H. is pure imaginary. Because, by the proposition, e^j 
is pure imaginary, it follows that Q2j is in the span of Ti and S. 



There are equations similar to (7.8) which give the evolution of the density gener- 
ating function, Q, under the higher flows. These have the form [FNR]: 

(7.11) Q tn = [Q2n-1,Q]. 

The loop Q has a natural decomposition Q = Q( ev ^ +Q( od ^ into, respectively, the even 
and odd terms of the series. For £ real, these are also, respectively, the imaginary and 
real components of Q. Moreover, these two parts are orthogonal with respect to the 
inner product, (■, •) on the loop algebra. This decomposition can be applied to (7.11) 
(for instance by equating imaginary and real parts on both sides of the equation). We 
apply this in the case of the odd flows; i.e., when n = 2m — 1 is odd. This yields two 
equations: 

Qtt 1 = [Q {ev \Q2 m -i], 
Ql:t 1 = [Q (od \Q2m-i] = 0- 

The right-hand side of the last equation is identically zero because [JF, JF] = 0. 

Thus the subspace of even densities is preserved under the odd flows. In particular 



(7.13) Q; 



1 2 1 a 
—u -eo x u 

2 2 

{ a 1 1 

-eo T u ——u 



which carries the densities that comprise the Miura transformation (2.7), v = u 2 + 
ied x u = 2(/i2 + e.2), preserves its form under the odd flows. In Section 2 we only 
considered the evolution of these densities under the £3-flow, the mKdV flow. 

Notice that for any of the elements in Q^ ev \ the coefficients of terms proportional 
to £ are formally O(e). Hence, in the dispersionless limit the odd NLS flows acting on 
real initial data collapse onto the Abelian subalgebra spanned by 7i which is comprised 
of the Hamiltonian densities introduced in Section 3. This extends to the higher odd 
flows our formal observations about the dispersionless limit of the mKdV flow which 
motivated the rigorous analysis of the zero-dispersion limit carried out in Sections 5 
and 6. 



8. Concluding Remarks 
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There is an alternative approach to showing that the limits of the focusing mKdV 
equation (1.6) are very close to those of the KdV equation (1.1). 

One is the Miura transform [Mr]. The simple relation 

— u = M(v) = v 2 ± ev x 

provides a link between the semiclassical limits of the KdV equation and of the defo- 
cusing mKdV equation: 

(u t + 6uu x + e 2 u xxx ) + (2v ± ed x )(v t - 6v 2 v x + e 2 v xxx ) = 0. 

This shows that the limit of v 2 , the mass density, (1.7), of the defocusing mKdV 
equation (1.6), converges to lim e ^ ' u ? the zero dispersion limit of the solution to the 
KdV equation (1.1), in the sense of distributions. This relation of semiclassical limits, 
between the defocusing NLS hierarchy and the KdV hierarchy, works at least in the 
case when their initial values depend only on the continuous spectrum of the associated 
ZS Dirac operator and the KdV Schrodinger operator Cs respectively [JLM2]. 

A similar transform, 
(8.1) u = -M(iv) = v 2 =F iev x 

gives a relation of solutions between a real valued focusing mKdV equation and a 
complex valued KdV equation: 

u t + 6uu x + e 2 u xxx = (2v =F ied x )(v t + Qv 2 v x + e 2 v xxx ). 

It suggests that the limit of v 2 , the mass density of the focusing mKdV equation , 

v t + 6v 2 v x + e 2 v xxx = 0, 

should converge to the zero dispersion limit of the KdV equation (1.1) as e — * in the 
distributional sense as we have in fact shown, for a reasonable class of initial data, in 
this paper. 
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u(x) for t=1.0 x-t-u contour plot 



Fig 2.1 Time-slices, u(x, t)|t=o.5,i.o,2.0j and contour plot (lower right) u(x,t) = const., 
< x < 27r,0 < t < 4. u(x, 0) = cos 2 x. The horizontal axis is x, and the vertical axis 
is either u (slice) , or t (contour plot) . 




Fig 4.1 The initial data A(x). Note its critical value r) max and the indicated defining 
relations for x ± = (77) . 




e = 0.02, t = 1.5 e = 0.02 contour plot 

Fig 2.2. Contour plot, u e (x,t) = const., < t < 1.5 (left). Slice, u e (x, t) |t=i.5, 
(right). Here, e = 0.08 (top), 0.04 (center), and 0.02 (bottom), and u e (x,0) = cos 2 x. 
The horizontal axis is x, and the vertical axis is either t (contour plot) or u (slice). The 
graphs have the same scales for comparison. 





d- 1 [<?M 2 - \uft 



Fig 2.3 Primitives of the first (left) and third (right) conserved densities. The cases for 
for e = 0.08 and e = 0.02, are overlaid, with the wider oscillations corresponding to large 
e. Here, the pseudo-inverse d" 1 is a periodic normalization of J x : each discrete Fourier 
mode of the density is multiplied by ^, except for mode 0, which is left unchanged. 
As e decreases, both the wavelength and the magnitude of the deviating oscillations 
decrease. This suggests that the limit for the primitives will be classical. 



